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Abstract 

Linked cluster expansions are generalized from an infinite to a finite 
volume on a d-dimensional hypercubic lattice. They are performed to 20th 
order in the expansion parameter to investigate the phase structure of 
scalar O(N) models for the cases of N = 1 and N = 4 in 3 dimensions. In 
particular we propose a new criterion to distinguish first from second order 
transitions via the volume dependence of response functions for couplings 
close to but not at the critical value. The criterion is applicable to Monte 
Carlo simulations as well. Here it is used to localize the tricritical line 
in a $ 4 + $ 6 theory. We indicate further applications to the electroweak 
transition. 

1 LINKED CLUSTER EXPANSIONS IN THE 
INFINITE VOLUME 



Convergent expansions such as linked cluster, hopping parameter or high tem- 
perature expansions provide an analytic alternative to Monte Carlo simulations. 
Originally they have been developed in the infinite volume, meanwhile they have 
been extended to finite volumes as well [El . In contrast to generic perturbation 
theory, hopping parameter expansions (HPEs) are convergent expansions about 
completely disordered lattice systems. The expansion parameter k is the co- 
efficient of the (pair) interaction term. Since we calculate free energies and 
connected correlations in the hopping parameter expansion, we generate linked 
cluster expansions (LCEs). 
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LCEs have been generalized to scalar field theories at finite temperature by 
. Compared to the zero temperature results, the finite temperature induces a 
tiny shift in the critical temperature (hopping parameter) . Thus one has to face 
a similar problem as the improvement of a bad signal to noise ratio in Monte 
Carlo simulations. Here the signal corresponds to the finite temperature effect 
and the noise to the "background" of temperature zero graphs. The price to 
describe the critical region, in particular critical exponents, with an accurracy 
of ks 1 % is to go to the 18th order in the expansion in k. This accurracy has 
been achieved in ||. Thus the question arises, why we are still interested in a 
generalization of LCEs to a finite volume. The reason is twofold. 

• The first one is to identify 1st order transitions from series expansions in 
the high temperature phase. For example a measurement of a vanishing critical 
exponent v ■ r\ is compatible with a Gaussian fixed point, but not conclusive 
for the onset of a first order transition region at a tricritical point. (The latter 
interpretation actually applies to our measurements.) 

• The second reason is to distinguish second order transitions associated 
with different universality classes. Ideally critical exponents should show a clear 
gap between plateaus of different universality classes. Practically the gap is 
smeared out because of the truncation of the series expansions at finite order. 
The truncation has a similar rounding effect on this gap as a finite volume on 
thermodynamic singularities. 



2 GENERALIZATION OF HPEs TO A FINITE 
VOLUME 

2.1 Embedding factors on the lattice 

o 

In HPEs the action is split into a sum of ultralocal parts S and a next neighbour 
part Snn with next neighbour couplings oc k. A Taylor expansion in k of the 
partition function Z about the ultralocal contribution to Z finally leads to 
graphical expansions of n-point susceptibilities \n 

The coefficients aj^' are a sum of graphs each of which consists of a product of 
the inverse symmetry factor, an internal symmetry factor, a lattice embedding 
factor, and a polynomial of vertex contributions depending on the couplings 

o 

involved in S- It is only the embedding factor that depends on the topology 
of the particular lattice, that will change in passing from an infinite to a finite 
volume. 
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2.2 Shift in the critical coupling 

From our series expansions we measure a shift in k c from a fit according to 

ln(|/c c (X) — k c (oo)\) ~ lnc — yx InL. (2) 

Here c is a constant, yx — 1/v for 2nd order transitions and yx oc d for 1st 
order transitions. The critical couplings k c (L) and k c (oo) are determined as 
the radius of convergence of the series expansions in a finite or infinite volume, 
respectively. Since an extrapolation of the scaling behaviour from the height 
and width of the critical region is not feasible within the series expansions, we 
have proposed the following monotony criterion. 



2.3 The monotony criterion 

Consider two volumes V\ and V% such that V\ < V2 < 00, < a < 1, a stands 
for a generic coupling, and t = (re — k c )/k c for the reduced "temperature". We 
define 



r(Vi,V a ):=l- Xa i; = a ?vV (3) 
X2(k = a ■ K C ,V 2 ) 



The monotony criterion says that 



> , 2nd order 
r Vi,v 2 { < 1 1st order 

= , tricritical point for dr/da ^ 0. 



The different behaviour of ry 1: y 2 comes from the singular part of \2- The 
underlying observation is a different monotony behaviour of response functions 
X for 1st and 2nd order transitions that is seen for %s with nonanalytic behaviour 
in the infinite volume limit. The %s are increasing in volume in a certain neigh- 
bourhood of T c for 2nd order transitions and decreasing for 1st order transitions. 
For a precise specification of the scaling region, i.e. the bounds that Vi, V2,\t\ 
have to satisfy, we refer to [|l] . The difference reflects the (5-function and power 
law singularities for 1st or 2nd order transitions, respectively, predicted at Tc 
in the thermodynamic limit. Thus the actual behaviour of xi we have found in 
applications to scalar O(N) models (cf. section 3) is neither a pecularity of the 
model nor an artifact of the series expansions. 

• It should be noticed that the criterion is not restricted to series expansions, 
but also applicable to Monte Carlo calculations (since both involved volumes 
may be finite) if the regular contribution to x is really negligible. 

• It is neither restricted to order parameter susceptibilities as the notation 
suggests, but similarly applies to other response functions like the specific heat 
if they diverge in the thermodynamic limit. 
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2.4 The effective potential for V < oo 



The effective potential is defined via the effective action T evaluated for a con- 
stant background field M (L d ■ V e ff := — T(M)). It can be expressed in terms 
of the lPI-susccptibilities according to 



1-ADk X \ pi M 2 xl PI (M 2 ^ 2 



Xl pi 2 ( X \ny 4! 



V* f(M) 



[x iPi f x ipi j 6! 

The X2« 7 s are obtained from the X2nS if the graphical expansion of the X2nS is 
restricted to IPI graphs. The (non)convex shape of the effective potential in 
the symmetric phase indicates the order of the transition. Its evaluation in a 
finite volume leads to an alternative determination of k c (L) via the coexistence 
of different minima of the free energy or a vanishing of the coefficient of the 
quadratic term. 



3 RESULTS FOR A3d $ 4 + $ 6 THEORY 

We apply the finite size scaling analysis with HPEs to a scalar theory with 4> 4 
and $ 6 -terms and O(N) symmetry in 3 dimensions, N is chosen as 1 and 4. 
After a suitable reparametrization the ultralocal part of the action reads 

S(<f>,\,a) = $ 2 + A($ 2 - l) 2 + cr($ 2 - l) 3 (4) 
The hopping parameter term is given as 

s nn = -2k ®( x ) ■ $ (y)- ( 5 ) 

x,y nn 

For a suitable choice of couplings the model has a tricritical line separating 
regions of first and second order transitions. Thus it serves as an appropriate 
toy model for testing our methods. 

A fit of the shift in k c according to EqJ| leads to the values of yr listed in 
the Table. 



N 


VT 




order 


1 


6.068(43) 


w 2d 


1st 


4 


5.55(59) 


w 2d 


1st 


1 


2.656(48) 


w 2jv 


2nd 


4 


not conclusive 




2nd 



The value for yx for the first order transitions confirms results of H for N = 1 
and extends results of pj for N = 4. It disagrees with the naive expectation of 
a scaling with d rather than 2d. The figure shows results for the monotony 
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criterion applied to the special case of Vi = 4 3 , V2 = oo 3 , a = 0.98, N = 4, 
A = er/2. The additional index M refers to the order of truncation. A similar 
figure is obtained for r(Vx, V2) as function of a for fixed order of the truncation 
and varying volumes. Thus the intersection of the curves with the cr-axis, which 
localizes the tricritical point, weakly depends on the order of the truncation 
and the choice of V\. The ultimate determination of the tricritical couplings 
at and At needs an extrapolation in both parameters. Respecting the mutual 
relation between the truncation effect and the finite volume dependence of the 
intersection, the final estimate for at leads to a t — 9.454(49) . For further 
details on the involved extrapolations we refer to This result improves 
the infinite volume estimate of 8 < at < 10 by two orders of magnitude. It 
should be noticed that the volume independence of \2 along the tricritical line 
supports the validity of a mean field analysis of tricritical exponents even in 
three dimensions as one might have expected from the Ginzburg criterion. The 
localization of at from the change in the shape of the effective potential leads to 
an upper bound on at, we find 9.75 < at < 10.0, since the involved systematic 
errors are less well under control. 

An extension to effective scalar models for an underlying Salam- Weinberg 
theory with HPEs in a finite volume will particularly well work in the range of 
large Higgs masses. It complements the range of Higgs masses which has been 
available in recent Monte Carlo simulations and allows a determination of the 
critical Higgs mass above which the electroweak transition turns into a smooth 
crossover phenomenon jj| . 
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